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Abstract. 

In most of the models of hadronic collisions the number of exchanged colour strings 
grows with energy and atomic numbers of the projectile and target. At high string densities 
interaction between them becomes important, which should melt them into the quark-gluon 
plasma state in the end. It is shown that under certain reasonable assumptions about the 
the string interaction, a phase transition to the quark gluon plasma indeed takes place in 
the system of many colour strings. It may be of the first or second order, depending on the 
particular mechanism of the interaction. The critical string density is about unity in both 
cases. In the latter case the percolation of strings occurs above the critical density. The 
critical density may have been already reached in central Pb-Pb collisions at 158 A GeV. 



US-FT/24-97 



*On leave of absence from the Department of High Energy Physics, University of St. Petersburg, 198904 
St. Petersburg, Russia. 



1 



2 



1 Introduction 

Multiparticle production at high energies is currently (and successfully) described in terms 
of colour strings stretched between the projectile and target [ 1-3]. Hadronizing these strings 
produce the observed hadrons. In the original version the strings live independently and 
the observed spectra are just the sum of individual string spectra. However, with growing 
energy and/or atomic number of colliding particles, the number of strings grows and one 
expects that interaction between them becomes essential. This problem acquires still more 
importance, considering the idea that at high enough energy collisions of heavy enough nuclei 
may produce the quark-gluon plasma (QGP). The interaction between strings then has to 
make the system evolve towards the QGP state. 

In our earlier publications we proposed a simple probabilistic model of string fusion [4], 
which was later realized as a Monte-Carlo algorithm [5]. It turned out to be quite success- 
ful in explaining some peculiar features of multiparticle production, such as heavy flavour 
enhancement at high energies and production of particles outside the kinematical limits of 
nucleon-nucleon collisions [6,7]. It did not however seem to predict any phase transition in 
the colour string system, which could be interpreted as a transition to the QGP. 

In a recent publication it was noted that the string density reached in Pb-Pb collisions at 
present energies is of the same order of magnitude as the critical density for the percolation of 
overlapping discs [8]. This could mean that strings may evolve to the QGP via the percolation 
process. It however remained not clear how the percolation can be incorporated into a more 
general picture of string interaction (and fusion), and why simple models of this interaction 
(in particular, the one proposed in [4]) do not lead to a phase transition. Also one would like 
to know possible alternatives to the percolation to be able to study from the experimental 
data the particular way the hadronic matter chooses to go into the QGP state. This is also 
related to the possible signatures of the percolation or of any other phase transition. 

This note intends to discuss these poroblems. Essentially our conclusions are as follows. 
We show that the model of Ref. [4] indeed lacks any phase transition. Its absence is due to a 
particular fusion mechanism chosen (Sec. 2). A slight change in this mechanism gives a very 
similar and also explicitly solvable model, in which there is a clear phase transition to the 
QGP. This transition is of the 1st order (Sec. 3). The percolation occupies an intermediate 
place between the two models in what concerns the interaction mechanism. Its predictions 
for observable quantities are discussed in Sec. 4. Sec. 5 contains some conclusions. 

This note is centered on the analytical treatment of the problem of string interaction. Such 
questions as string formation and their hadronization, which require Monte Carlo simulations, 
remained outside its scope. As a consequence we shall study only very primitive observables, 
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whose behaviour under string fusion can be predicted on plausible assumptions: the total 
(charged) multiplicity and the total number of strings. The latter shows itself in forward- 
backward correlations [9]. We shall study the behaviour of the system of a fixed number N of 
strings under their interactions. In practice, of course, the number N fluctuates, its average 
value growing with energy. This will undoubtedly smooth all irregularities introduced by the 
phase transition and make its detection still more difficult. 



2 Fusion of strings in the model [4] 

In this section we briefly review the model proposed in [4] from a somewhat novel point of 
view, looking for a possible sign of a phase transition in it. It was assumed in the model 
that strings fuse as soon as their positions in the transverse space coincide. The probability 
of fusion p was assumed to be the ratio of the string own transverse area v to the total 
interaction transverse area V. Of course one should consider v as an effective area, which 
may be smaller than the geometrical one for small fusion probability. We also assumed that 
colours of fusing strings (i.e. colours placed at their ends) sum into the total colour of the 
resulting string. For the realistic QCD this summation should be performed according to the 
composition laws of the SU (3) colour group. This can only be done within a Monte Carlo 
code. In the analytical approach of [4], as well as in the present note, to see qualitative 
predictions of string interaction, we adopt a simplified, Abelian, approach, taking that two 
strings with positive integer colours n\ and ri2 fuse into a string of colour ri\ + n<i- Ordinary 
string possess a unit colour n = 1. The basis of the model is laid by the probability P(v n ) to 
find v n strings of colour n, chosen to be 

/~\N N-M M—l 

Here N = J2 n nv n is the fixed total colour of the strings which is also the number of ordinary 
strings originally produced; M = J2 n u n i s the total number of strings ("clusters", in terms 
of the statistical approach). In the following we pretend to change this probability. So it is 
instructive to analyse its foundation. 

Since we shall study the behaviour of the system at fixed N, factors depending only on N 
do not play any role and are actually a part of the normalization constant. The denominator 
factors in (1) are of the standard quantum statistical origin: v n \ takes into account the 
identity of v n strings of colour n and each n! takes into account the identity of n ordinary 
strings which fuse into a string of colour n. The number N — M is, in fact, the number of 
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fusions, hence the factor p N . We come finally to the factor 

M-l 

w= H(l-kp) (2) 

k=l 

which is the one that contains information about a particular way strings fuse. 

To obtain it we made an additional assumption that fusion does not change the proper 
volume v of strings, that is, all strings have the same volume v irrespective of their colour 
(and consequently the "string tension" of a string with colour n is n 2 greater than of the 
ordinary string). As a result, fused strings have the same probability to fuse as ordinary 
ones. The factor w can then be obtained from a simple recurrency relation, which follows if 
one studies the system with one string of arbitrary color k added. For this system to exist, 
it is necessary that the new string should not fuse with the rest. This leads to a relation 

w{v n + 5 nk ) =w(v n )(V-vJ2 u n)/V (3) 

n 

The second factor, equal to 1 — pJ2 n v n, gives precisely the probability not to fuse with the 
rest. We chose to write it in terms of volumes to compare with the treatment in the following 
section. The recurrency relation (3) has evidently the solution (2) 

Given the probability (1), one can calculate all observables of interest. Our main observ- 
able will be the average number of strings of colour n, (u n ). One easily finds (see Appendix 
for some details) 

K> = E V nP(yn) = C n N p n -\l-p) N - n (4) 

It is just a binomial distribution in n. With the help of (4) we obtain the average total 
number of strings (M) : 

(M) = ±(l-(l-p) N ) (5) 

From this last equation we see that (M) monotonously falls from (M) = N at p = to 
(M) = 1 at p = 1, which clearly shows the effect of the fusion process but exhibits no sign of 
irregularity, which could indicate a phase transition. 

One would also like to know the behaviour of the average multiplicity (/i) of particles 
produced by the hadronizing strings as a function of p. To this purpose one has to make 
certain assumptions about the multiplicity fi n corresponding to a string of colour n. On 
physical grounds the multiplicity should be proportional to the surface of the string area. 
With all strings having the same dimension this does not lead to any n dependence. One 
may also expect some growth of the multiplicity with n from the growing string tension. In 
the limiting case when the string of colour n emits particles exactly as n ordinary strings we 
shall have 

n 
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that is, one will not see any change in the spectra resulting from fusion. This is not a surprise, 
since with p n = n\i\ fusion does not influence the multiplicity. In the opposite case, when a 
high-colour string emits particles exactly like an ordinary one, p n = p\ and one has 

(fi) = mM 

and with the growth of p (fi) falls from Nfii to p\. 

More generally, one can study the distribution in the number m of produced particles, 
assuming that this distribution is A n (m) for a string of colour n, satisfying 

A n (m) = 1; ^2 m ^n(m) = p, n 

m m 

The resulting distribution A(m) is then given by the average 

A(m) = E ^ H < Vn > (6) 

En < Vn > 

From (6) one can calculate moments of the distribution A(m). In particular the average 
number of emitted particles per string is 

< m >=< fi > I < M > (7) 

and the dispersion squared is 

D 2 = Y2 m 2 A(m)- < m > 2 (8) 

m 

Presenting 

Y,m 2 \ n {m) =n 2 n + d 2 n (9) 

m 

where d n is evidently the dispersion of the distribution in m for a string of colour n, we 
split D 2 into an internal part, generated by this intrinsic dispersion, and an external one, 
generated by fluctuations in the distribution of strings in colour. This latter is given by 

Dl x = < Vn > ~ < m > 2 (10) 

n 

and is presumably sensitive to a possible phase transition. It is positive and vanishes both 
for p = and p = 1. 

To have a clearer vision of the model from the statistical point of view, it is instructive 
to consider the "thermodinamical limit", which in our case means V — > oo and N — > oo with 
the density of strings p = N/V constant. Evidently then p = v/V — > and the relevant 
parameter describing string interaction is 



rj = pv = Np 



(11) 
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In this limit we find from (4) and (5) 



{Vn)/N = 



1 



exp n — r) + n In — 



) 



(12) 




and 



(M)/N = - (1 - ei) 



(13) 



We observe that the model has correct scaling properties with N: both {v n )/N and (M)/N 
result independent of ./V in the limit. From (12) we conclude that (u n ) has a maximum at 
n = r], so that with the growth of rj the distribution steadily shifts towards higher n. However 
there is no phase transition: the system does not change qualitatively as i] — > oo. 

This circumstance is clearly seen from Figs. 1 to 3, which show < M >, < fi > (in units 
fii) and D\ x (multiplied by 10) as a function of the parameter rj in the system of 50 ordinary 
strings. In this figures prediction of the model studied in this section are illustrated by the 
upper curves (marked 1). The curves are quite smooth and show a very slow fall off of < M > 
and < /i > and an equally slow growth of D 2 ex with the growth of r]. 

The reason of this is, of course, evident. It resides in the assumption that the dimension 
of strings does not grow with colour. Strings with higher colour then fuse with the same 
probability as ordinary ones, so that fusion is not accelerated with the accumulation of high 
colour strings. This picture does not correspond to what one might expect on intuitive 
grounds. Presently we shall see, that once we substitute this assumption by a more physical 
one, the behaviour of the system changes radically and acquires the expected features with 
a clear phase transition at high enough p. 

3 Fusion of strings with the conservation of their volume: 1st 
order phase transition 

An alternative model, which also preserves the important quality of being explicitly soluble, 
is obtained if we assume that as string fuse, they conserve their proper volume. Namely, two 
ordinary strings fuse into a string with the volume 2v. Then generally a string of colour n 
will have volume nv, and consequently the probability to fuse much higher. To obtain the 
relevant factor w we construct a recurrency relation similar to (3). However in this case the 
factor w will depend not only on the occupation numbers v n but also on the volume V, since 
the probability of fusion depends on the volume available and the latter now depends on the 
colour of the added string. The recurrency relation we obtain is 



w 



(V 7 , u n + 5 nk ) = cw(V - kv, v n )(V -kv-v^Vn) 



(14) 



n 
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We have taken into account that the rest strings live in a smaller volume V — kv left after a 
string of colour k has been added. 

Recurrency relation (14) can also be solved quite easily. We rewrite it as 

w(V + kv,u n + S nk ) = cw(V,v n )(V -v^v n ) (15) 

n 

Note that the factor V — v J2 n u n is invariant under a simulataneous change V — ► V + kv, 
v n — ► + ^nfe- Therefore, making such substitutions we obtain 

w(V + fc gfc v, v n + % 5 nfe ) = c 9fe u;(F, u n )(V - v ^ ^) 9fc (16) 

and more generally 

^ + «E fc 9*' 1/ '» + ft») = cE 'N^''n)^-1'E^ E * ft ( 17 ) 
Putting finally v n = 0, io(V, 0) = 1 and c = 1/V we obtain 

w = (l-p7V) M (18) 
Thus the distribution in v n which desribes our new model is 

n N n N-M 

One observes at once that this distribution makes sense only for p < 1/N, or, in terms of 
the parameter i], Eq. (11) 

V < 1 (20) 

What happens if p becomes greater than 1/7V? The answer is evident. The total volume 
of strings Nv then exceeds the interaction volume V and the only way for the strings to 
exist is to form a single string which spans all the interaction area. In other words, with the 
probability one, all string fuse into what may be viewed as the QGP plasma. This transition 
is of the first order, as we shall presently see considering the thermodynamic limit. 

With the distribution (19) one can also calculate average occupation numbers in the 
(nearly) explicit form. One finds (see Appendix) 

K> = NaC- N l^§^ (21) 



and 



where 

a = I/77 — 1 (23) 
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and /at (a) is denned to be 

(24) 

x=0 

It is a polinomial of order N. 

These exact expressions are not very transparent. However one clearly sees from (22) 
that with the growth of rj from to 1 (that is p growing from to 1/N) the average value 
(M) diminishes from N down to exactly 1. Indeed, at rj — > we have a — > oo and only two 
highest powers of the polinomials /n have to be taken into account. One finds 

f N (a) ~ a N + (1/2)N(N - l)^' 1 + ... 

Putting this into (22) one obtains (M) = N. At rj — > 1 we have a — > and only the lowest 
term of f n contributes, which is a independently of N. From (22) one then finds (M) = 1. 

To clearly see the phase transition we again take the thermodynamic limit V, N —> oo with 
the density p and parameter rj fixed. Applying the saddle point technique (see Appendix) we 
then find a simple formula for (u n ) 

{v n )/N = au n /n\ (25) 

where a is given by (23) and u is a solution of the transcendental equation 

au = e-' a (26) 

Eq.(25) is nothing but a Poisson distribution (correctly normalized because of (26)). For the 
average (M) we obtain 

(M)/N = l/u-a (27) 

As we see, the system, again, possesses normal scaling properties with respect to the volume: 
the right-hand sides of (25) and (27) are independent of N. 

From (25) one concludes that there is no infinite cluster for rj smaller than the critical 
value rj c = 1, since = 0. At rj = rj c , a goes to zero and u becomes infinite, so that Eq.(25) 
ceases to exist. At rj > rj c we have a single infinite cluster and no finite ones. Thus if Pf (oo) 
are the probabilities to find a finite (infinite) cluster we have 

P f = l, Poo = 0, < v < Vc = 1 

P f = 0, Poo = 1, r] > Vc = 1 

This situation corresponds to a 1st order phase transition at rf = rf c 
phase transition point we have 

(M)/N = 1 a^i-n^O 

| In (aln-1 | 



/at (a) 



d_ 

"dx' 



-i) 



1. Slightly below the 



(28) 
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so that the behaviour near the critical point is logarithmic, and not power-like, as for the 2nd 
order phase transitions. 

As to the multiplicity, it is reasonable to assume that it is proportional to the emitting 
surface and therefore to \Jn for the string of colour n. Then the observed average multiplicity 
will be 

<M>=Mi£v^<"n> (29) 

n 

Hi being as before, the multiplicity for a single ordinary string. The sum (29) cannot be done 
analytically for finite TV. For large N the saddle point method gives 

{fj)/N = W /v^ (30) 

where u is the solution of (26). This expression is valid for not very small u: evidently for 
rj = and consequently u = we have 

(fi)/N = W (31) 

At rj close to the critical point we have 

(H)/N = , ^ . (32) 



In (a In ±) 



logarithmically small, as the number of strings. Above the phase transition the whole inter- 
action area acts as a single emitter and we have 

M/N = (33) 



so that the multiplicity is down by a (large) factor y/~N. 

These features of < M > and < /x > are clearly seen in Figs. 1 and 2 (curves marked 
2). One finds a sharp drop of both < M > and < /j, > at the phase transition point 77 = 1. 
Substantially below the transition point the fall off of both these quantities is quite slow 
and for 77 below 0.5 predictions of the model practically coincide with those of the preceding 
section. The same conclusions can be drawn from Fig. 3 where the curve marked 2 shows 
the dispersion squared D\ x . Its sharp rise as 77 approaches unity and subsequent fall to zero 
at 77 exactly equal to unity are quite indicative of the first order phase transition. 



4 Percolation 

An alternative way to the formation of the QGP in the colour string system is the percolation. 
This is a purely classical mechanism. Consider projections of ordinary strings onto the 
transverse space as filled circles (discs) of area v. Let them be distributed uniformly in the 



10 



interaction area V. Some of the discs may overlap and in this way form clusters of string 
matter. At some critical disc density p a cluster appears which spans all the volume (an 
infinite cluster if the number of discs N — ► oo). This infinite string cluster may be identified 
with the GQP filling part of the total volume. 

Comparing this mechanism with the string fusion we note that it occupies an intermediate 
position between the two scenarios discussed in Sees. 2 and 3, as far as the interaction 
mechanism is concerned. Indeed, in the percolation, strings are allowed to partially overlap, 
so that the volume of the resulting clusters varies between the two extremes: the volume of 
a single string, in the complete overlap case, and the sum of the interacting strings volumes, 
in the minimal overlap case. The two scenarios discussed above correspond just to these two 
extreme cases. An additional advantage of the percolation treatment is that it introduces 
various clusters for a given high colour n > 1, thus taking into account (in a classical, not 
quantum way) possible excited states of fused strings, which should become numerous for 
large n. Unfortunately, the percolation cannot be treated analytically and nearly all the 
results known in the literature have been obtained by Monte Carlo simulations in the vicinity 
of the phase transition point. 

The percolation phase transition is known to be of the 2nd order. The critical value of 
the parameter r\ introduced above is found to be [10] 

ilc ^ 1.12 (34) 

This value is remarkably close to the critical value r\ c = 1 for the 1st order phase transition 
in our model of string fusion discussed in Sec. 3. It is also known that the fraction <f> of the 
volume, occupied by strings is determined by the formula [11] 

4>=l-e- T > (35) 

which coincides with the fraction of the volume <fi = v(M) in our string fusion model of Sec. 
2 (Eq. (13)). In the model of Sec. 3, of course, 4> = 1, since the volume is conserved in this 
model and, in terms of percolating discs, no overlapping takes place. 

Below the phase transition point, for rj < 7] c , there is no infinite cluster, so that 

Pf = (1/A0 J2 niJ n = 1, ^oo = (36) 

n 

Above the transition point, at rj > r] c an infinite cluster appears with a probability 

Poo = 0( V - Tfc) (f? - 7] c f (37) 

The critical exponent (3 can be calculated from Monte Carlo simulations. However the uni- 
versality of critical behaviour, that is its independence of the percolating substrate, allows to 
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borrow its value from lattice percolation, where [10] 

/3 = 5/36 

The main quantity of interest to us will, as before, be the average occupation numbers 
(v n ). Their behaviour at all values of ry and n is not known. Prom scaling considerations in 
the vicinity of the phase transition it has been found [12] 

(u n ) = n- T F(n a (r) - r/ c )), \ V - Vc \ « 1, n » 1 (38) 

where r = 187/91 and a = 36/91 and the function F(z) is finite at z = and falls off 
exponentially for \z\ — ► oo. Expression (38) is of limited value, since near rj = r\ c the bulk 
of the contribution is still supplied by low values of n, for which (38) is not valid. However 
from (38) one can find non-analytic parts of various observables at the transition point. In 
particular, we find the singular part of (M) in the form 

A(M) = c\q - Vc\ 8/3 (39) 

This singularity is quite weak: not only (M) itself but also its two first derivatives in r\ stay 
continuous at r\ = r/ c and only the third blows up as \rj — r/d^ 1 / 3 . 

A little stronger singularity is found in the average multiplicity. Assuming it again to be 
proportional to linear dimensions of strings in the transverse plane, that is, to \pri, we get 
from (38) 

A( /U ) = C |t ? -7 ?c | 101 / 72 (40) 

So already the second derivative of (fi) blows up at the transition point. 

More detailed predictions for the total (M) and (/i) , and not only for their singular 
parts, can only be obtained from Monte Carlo simulations. For a system of 50 strings they 
are presented in Figs. 1 to 3. by the curves marked 3. One observes that the behaviour of 
both < M > and < fi > is quite smooth and does not seem to show any sign of irregularity at 
the phase transition point rj ~ 1.12. This is of no wonder, considering very weak singularities 
(39) and (40) introduced by the phase transition. Nevertheless, on the whole, the behaviour of 
both quantities is very different comparing to the previous models. Below the phase transition 
the string fusion is much stronger in the percolation case, so that averages of both M and /i 
result much smaller than in the fusion models of Sees. 2 and 3. Above the phase transition 
point these averages continue to fall off quite smoothly, so that they become substantially 
greater than for the 1st order transition (Sec. 3) but stay much smaller than in the model of 
Sec. 2 without a phase transition. This opens up some possibilities to study the interaction 
mechanism from the observed behaviour of < M > and < fi >. The dispersion curve in Fig. 
3 is more informative. It shows a clear maximum at values of r] in the interval 1.1-1.5. The 
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absolute value of the dispersion at the maximum is nearly two orders of magnitude higher 
than those for the fusion models of Sees. 2 and 3 (which are shown in Fig. 3 multiplied by 10). 
Such a high dispersion is typical for a phase transition. From Fig. 3 one can approximately 
determine the critical value rj c to lie around 1.3. It is somewehat higher than the value 1.12 
known in the literature, but the latter corresponds to infinitely many strings, and our results 
refer to only 50. 

5 Conclusions 

We have shown that string fusion may lead to a phase transition into the QGP for high enough 
string density. We presented a particular completely solvable string fusion model in which 
this phase transition is of the 1st order. Percolation is another alternative, which has certain 
advantages in treating the details of string interaction. It is important that in both models 
the transition to QGP starts practically at the same value of the string density (corresponding 
to rj = 1 in the string fusion model and r] = 1.12 in the percolation model). So, irrespective 
of the interaction mechanisms, one can expect QGP formation at such densities. As pointed 
out in [8], this critical density seems to have been nearly reached in Pb-Pb collisions already 
at present energies at SPS and should be considerably surpassed in future RHIC and LHC 
experiments. 

We have dealt only with qualitative features of string interactions and transition to QGP, 
which allow for an analytical treatment. Detailed predictions require a Monte Carlo algo- 
rithm. 

The two models of phase transition presented above seem to be complementary for such a 
program. The string fusion model of Sec. 3 is well suited for the hadronization part, since one 
can comparatively easily extend the standard formalism of string decay to strings of higher 
colour. In fact, this already has been included in the Monte Carlo string fusion scheme of Ref 
[5]. However the model is difficult to realize on the level of string production: somehow the 
original strings should be melted into the fused one of the same geometrical form, which is 
not at all easy in the essentially classical Monte Carlo approach. The percolation treatment, 
on the contrary, is quite easily and naturally realizable at the string formation stage, but 
meets with serious problems at the hadronization level, when one has to introduce the decay 
of complicated geometrical structures corresponding to percolating clusters. These problems 
are left for future studies. In any case, experimental signatures like J/* suppression in Pb-Pb 
collisions, seem to indicate that percolation of strings already occurs [8]. 
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8 Appendix: Averages over the distributions (1) and (14) 

We illustrate our technique by calculating the normalization constant c in the distibutions 
(1) and (19). We use the identity 

dz 



-Sir/" < 41 » 



where the integration is supposed to go along a closed contour around the origin. Then for 
both distributions we obtain 

du _ M f dz ^ uEn^^En"" 



c- 1 = (QP) N £ / tF-«- N I —z- M P- M w(N, M) V JT" "f" " (42) 

where the sums over v n are now unconstrained and can be easily done. In fact, we have 

Tn^^-^=eM^"/n» (43) 

and the product over n gives 

T=l[T n = ex P (z(e u -l)) (44) 

71=1 

Integration over z is trivially performed and gives 

(e u - l) M /M\ 

Now we have to sum over all M. This summation is different for the two distributions 
(1) and (19). For (1) the factor w is given by (2) and can be presented in the form 

W {M) = v M M\Cf jv 

Then the sum over M gives 

e u ' p (45) 

For (19) the factor w is given by (18). The summation over M then gives a more complicated 
factor 

e iV«(e«-l) (46) 

where a is given by (23). 

Finally the integration over u has to be done, which reduces to taking the iVth term of 
the expansion of expressions (45) or (46) around u = 0. For the distribution (1) it is trivial 
and gives 

c- 1 = Q N /N\ (47) 
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whereas for the distribution (19) one obtains 

c- 1 = (Qp) N f N (Na)/N\ (48) 

where Jn is defined by (24). 

Note that the polinomials /at (a) can be easily calculated by a recurrency relation which 
relates the coefficients for and /jv+i- Indeed one notices that the derivatives of the function 
(46) have the structure 

{d/du) N e a{eU ^ = e a{eU ^f N (ae u ) (49) 
Differentiating once more, one gets a relation 

cP +1 =c%- 1] + kcP,l<k<N (50) 

for the coefficients cffl of the polinomial fx. Also one evidently has c$ = and cffl = 1. 
Thus, starting from N = 1, one can successively calculate all polinomials /at using the 
recurrency relation (50). 

Averages like (v n ) or (M) have the same representation as (42) with some simple addi- 
tional factors and are calculated in a similar manner. 

To study the limit N —>■ oo the integral representation for /jv(ATa) is more convenient: 

= / 2^»-^" (e "- 1) < 51 » 

Presenting the integrand in the form exp NE{u) we have at large N 

E~a{e u -l)-lnu 

so that the saddle point is determined by Eq. (26). Knowing the saddle point it is trivial 
to calculate the averages {f n ), since for them an extra factor Nau n /n\ has to be introduced 
under the integral. Taking this factor at the stationary point one obtains Eq. (25). 
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9 Figure captions 

Fig.l. The average number of clusters M in the system of 50 strings as a function of 77. 
Curves 1, 2 and 3 correspond to models of Sees. 2, 3 and 4 respectively. The latter curve 
was obtained by Monte Carlo simulations. 

Fig.2. The average multiplicity /x in units \i\ in the system of 50 strings as a function of 
rj. Curves 1, 2 and 3 correspond to models of Sees. 2, 3 and 4 respectively. In all cases fi n 
was assumed proportional to y/n. The percolation curve (3) was obtained by Monte Carlo 
simulations. 

Fig. 3. The dispersion squared D\ x in the system of 50 strings as a function of r/. Curves 1, 
2 and 3 correspond to models of Sees. 2, 3 and 4 respectively. Curves 1 and 2 show values 
of D\ x multiplied by 10. Curve (3) was obtained from Monte Carlo simulations for < M > 
and (fi) (Figs. 1 and 2). 




Fig. 1 



